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(X) \ We study the nine possible rare charm meson decays D Pl^l~ (P = vr, rj, rj') using 

^ ■ the Heavy Meson Chiral Lagrangians and find them to be dominated by the long distance 

^ . contributions. The decay tt^I^I^ with the branching ratio ~ 1 x 10~^ is expected to 

1^ ! have the best chances for an early experimental discovery. The short distance contribution 

in the five Cabibbo suppressed channels arises via the c — ul^l~ transition; we find that this 

\ contribution is detectable only in the D — >■ 7rl~^l~ decay, where it dominates the differential 

^ ■ spectrum at high-g^. The general Minimal Supersymmetric Standard Model can enhance 

T— I ■ the c —>■ ul^l~ rate by up to an order of magnitude; its effect on the D Pl^l~ rates is 

^ . small since the c ul^l~ enhancement is sizable in low-g^ region, which is inhibited in the 

' hadronic decay. 
Ph. 

ix. 

^ ■ 

■ 1 Introduction 

I The flavour-changing neutral processes are rare in the standard model and are of obvious 
interest in the search for new physics. Processes like c ^ wf and c ul^l~ are screened 
by the long distance contributions in the decays of charm hadrons JI], 0] and one has to look 
for speciflc hadronic observables |]^, |, Q in order to probe possible new physics [7|. The 
long distance contributions are expected to dominate over the short distance contributions 
also in the — mixing p|, for which interesting experimental results have been reported 
recently 

The long distance and the short distance contributions to rare charm meson decays 
D Vl^l~ with V = p,uj,(f), K* have been considered in 0. The long distance contributions 
were shown to be largely dominant and screen possible effects of new physics in c — > ul~^l~ , 
unless these are very large. The experimental upper bounds on their branching ratios are 
presently in the 10~^ range and are an order of magnitude larger than the standard model 
prediction for speciflc channels 0. The decay D'^ — >• p^l^l~ is predicted at the highest rate 
~ 3 X 10~^ P], but there is unfortunately no experimental data on this particular channel. 
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In the present paper we consider the weak decays D Pl^l^ with pseudoscalar P = 
TT, K, 7], 7]', some of which having contribution from the c — > ul^l~ transition. These channels 
have not been observed so far and only experimental upper bounds on the various branching 
ratios in the range 10~^ — 10~^ exist [0, |T2|, |l3l. The recent E791 analysis considers all 



and decay channels. The most recent FOCUS analysis [|T^] provides upper bounds 
of about 8 X 10~^ on the 7i^fi^fi~ and K^fi^fi~ branching ratios and is not 

far from our standard model prediction 1 x 10^^ for n^fi^fi^ . The limits on and 

modes at the level 10~^ are expected from CLEO-c and B-factories, while the limits on 



modes are expected to be an order of magnitude milder |n 



On the theoretical side, the long distance contributions to ^ ivl^l decays have been 



considered in |T^. We consider here also the long distance weak annihilation contribution 
and confirm it to be small in this channel. Calculations for other D Pl^l~ channels are 
not available in the literature. In the present work we investigate all these channels, including 
long-distance (LD) and possible short- distance (SD) contributions arising from the c — >■ ul^l~ 
transition. The QCD corrections to c — > ul^l~ amplitude have not been studied in detail 
yet and we incorporate only what we believe to be the most important QCD effects. We 
explore also the sensitivity of c — > ul^l~ transition to (i) minimal supersymmetric model 
with general soft-breaking terms and (ii) two Higgs doublet model with flavour changing 
neutral Higgs interactions. 

The c — > ul^l~ transition in SM, MSSM and Two Higgs Doublet model is studied in 
Section 2. The long distance contributions are considered within the Heavy Meson Chiral 
Lagrangian approach in Section 3. The results are compiled in Section 4, while conclusions 
are given in Section 5. 



2 The c ul^l decay 



The Lagrangian leading to c — ul^l transition is (using notation as in |16 



where 



-m^ua'"'PRC F. 



167r2 



O'^ = UJ^PrC h^l O[o = U-i^PRC h^j^l 



(1) 



(2) 



with Pji L = (1 ±75)72. In Eq. (|T|) only the CKM matrix element V^*Kis appears, for reasons 
explained in the next subsection |T^. The Wilson coefficients in various scenarios are given 
in the following sections. The differential branching ratio is given by |T6| 



dBr{c ul+l-) 
ds 



1 dT{c ul+t 



ds 



1927r3r(DC 



a 

4^ 



\v:,Vus\\i-sf 



{(1 + 2s)(|c9p + IcioH + 4(1 + 2/s)|c7p + l2Re[c*c^]] + {07,9,10 ^ ^7,9,1^] 



(3) 
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where s = mfi/ml, rric — 1.5 GeV and the mass of / = e,/i is neglected. The short-distance 
part of the D Pl^l~ amphtude, which is induced by c ^ ul^l~ transition, is given by 
(p3|) in Appendix. 



2.1 Standard model 



The c ul^l^ amphtude is given by the 7 and Z penguin diagrams and W box diagram at 
one-loop electroweak order in the standard model, and is dominated by the light quarks in 
the loop. One has @, |17 



4 nis 
- In — 

9 WLd 



1.34 ±0.09, C7,io(mvy) oc 



m 



d,s 







m 



w 



4,9,10 ("^14^) OC — C7,9,io ^ 



(4) 



for irLs/md = 21 ±4 MeV |]T3[, where the terms proportional to s/''^w have been neglected. 
The leading term ln(ms/mrf) in cg arises from the penguin diagram with photon emitted from 
the intermediate quark. 

The QCD corrections to c — ul^l^ amplitude have not been studied in detail yet. The 
QCD corrections to C7, which is extremely small at the one-loop level, have been studied in 
m and are found to be large 



nir. 



-(0.007 + 0.020 i)[l±0.2] . 



(5) 



We expect the QCD corrections to cg to be rather unimportant, given that cg is relatively 
large already at one-loop level [|l9l. We assume that the QCD corrections to cio do not 
affect the c — > ul^l~ rate significantly and use therefore only the C7 and cg coefficients. The 
differential branching ratios for the cases with and without QCD corrections are shown by 
solid and dashed lines in Fig. respectively. The branching ratio [6± 1] x 10"^ is small and 
arises mainly from cg; the contribution from C7 is small in spite of QCD enhancement. 

dBr ds 




0.2 0.4 0.6 



0.8 



Figure 1: The differential branching ratio dBr{c —>■ ul^l~) / ds: the dashed line denotes 
the one-loop standard model prediction, while the solid line incorporates also the QCD 
corrections to C7 [|l^. The best enhancement of the c — >■ ul^l~ rate in the general MSSM is 
given by the dot-dashed line, where the mass insertions are taken at their maximal values 
(I, 0) and as = 0.12, M^q = Mgi = 250 GeV. 
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2.2 Minimal supersymmetric standard model 

New sources of flavour violation are present in the minimal supersymmetric standard model 
(MSSM) and these depend crucially on the mechanism of the supersymmetry breaking. 
The schemes with flavour-universal soft-breaking terms lead to contributions proportional 
to g jj V*q^uq^l and have negligible effect on the c ul^l~ rate Our purpose here 



is to explore the largest possible enhancement of the c ul^l" rate in general MSSM with 
non-universal soft breaking terms. Based on the experience from the c — > U7 decay 
0, where the dominant contribution arises from gluino diagrams with the squark-mass 
insertion (^12)^^5 concentrate only on the gluino exchange diagrams with single mass 
insertion|^. Following the analogous calculation for b sl^l~ |T^, we get for the Willson 
coefficients in the MSSM 

4™"° ^ (8) 

with as = asimw) = 0.12, N, = 3, z = M^JM^^, Pi,k{z) = dx dy y'il-yy[l-y + zxy + 
z{l —x)y]~'^, e„ = 2/3 and = —1/3. The numerical bounds in (|^,0) are obtained by using 
parameter values as discussed below. The expressions for C7 9 are obtained by replacing 
L ^ R in the formulas above. We use gluino mass Mgi = 250 GeV and the common value 



for squark masses Msq = 250 GeV, given by the lower experimental bounds p3 |. 

The mass insertions are free parameters in a general MSSM. The strongest upper bound 
on {S'^2)lr comes by requiring that the minima of the scalar potential do not break charge 
or color, and that they are bounded from bellow |22, ^, giving 



KIlr, \S'^2\rl<0.00A6 for = 250 GeV. (9) 

The insertions (5i2)ll and {S12) rr can be bounded by saturating the experimental upper 
bound Am/5 < 4.5 x 10~^^ GeV by the gluino exchange |2^; the corresponding con- 
straint on {Si2)lr is weaker than @. Since we are interested in exhibiting the largest possible 
enhancement of the c ul~^l^ rate, we saturate Amo by {S^2)ll, obtaining [|], ^ 

\5i2\ll <0m for = Mg; = 250 GeV (10) 

and set {5i2)rr = 0. 

The biggest possible enhancement of the c — > ul^l~ rate is obtained using the mass 
insertions at their upper bounds and is shown by the dot-dashed line in Fig. The effect 
is dominated by the gluino exchange diagrams induced by {812) lr. and can enhance the 
c — >■ ul^l~ rate by nearly an order of magnitude, with the best enhancement displayed in 
Table |1|. 

^We work in the super-CKM basis for squarks, where the squark - quark - gaugino vertex has the same 
flavour structure as the quark - quark - gauge boson vertex; for review see pl[ . 
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The supersymmetric enhancement of c — > ul^l~ is due to the increase in C7 (Eq. ^ and 
is manifested at small mu due to the exchange of an almost real photon. This enhancing 
mechanism is unfortunately not present m. D ^ Pl^l~ decays (see Eq. (|23|)) since the decay 
D — > P7 with the real photon in the final state is forbidden (see Eq. (|13D). 







best enhanc. 


c ue^e 

c UjJ,~^l^~ 


(6± 1) X 10-9 
(6± 1) X 10-9 


6 X 10-s 
2 X 10-s 



Table 1: The second column represents the standard model prediction for c — *■ ul^l~ 
branching ratios, which is practically unaffected by the QCD corrections (see text). The 
third column represents the biggest possible enhancement of the branching ratio in MSSM, 
evaluated for mass insertions at their maximal values (|^, p!0|) . 



2.3 Flavour changing neutral Higgs 



The tree-level exchange of flavour changing neutral Higgs turns out to have a negligible 
effect on c — > ul^l~ rate, due to the strong constraint coming from the experimental upper 
bound on Am/) and due to the small mass of the leptons e and /x. Assuming the same 
c — u — H coupling^] f^u and mass mu = 300 GeV for all three neutral physical Higgses in the 
Two Higgs Doublet Model, and saturating the experimental upper bound Am^i < 4.5 x 10"^^ 
GeV m 



4/, 



2 

CM p2 



we get feu < 2 X 10 . This leads to a branching ratio 

Br{c - uf^^f^-r = ;;777^SLv ( ^) ' < 7 x 10-^^ . (12) 



7687r3r(D0) V vmjj 

Thus, unlike in the supersymmetric model, the experimental upper bound on ArriD imposes 
this new contribution to be negligible. 



The authors of have studied the constraints on the parameters of this model imposed 
by the present data on the semileptonic and leptonic D decays. Since they did not consider 
the constraint coming from the — mixing, they have obtained rather mild constraints. 



3 Long distance contributions 

Now we turn to an estimate of the long distance contributions to the D —>■ Pl^V decays. 
The dominant long distance contributions arise via the weak transition D — > P7* followed 

^The coupling is feu for c — u — ^ and fcul5 for c — u — A*^. 

■^The matrix elements of four-fermion operators is evaluated according to p^. 
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by 7* — > l^l . The general Lorentz structure of the D P'j* amphtude, consistent with 
electro-magnetic gauge invariance, is 



A[D{p) ^ P{p'h*{q,e)] cc A{q') [q\p + pT ' « - mi)q^] (13) 

and this amplitude vanishes for the case of a real photon. The factor q^ in (|1^) cancels the 
photon propagator and the general amplitude has the form 

A[D{p) ^ P7* ^ P/+(p+)/-(p_)] =i%e' A{q') u{p.)MP+)- (14) 

The long distance contribution is induced by the effective nonleptonic weak Lagrangian 

^|Ac|=i ^ [ai W7^(l - 75)g. qMl - 75)c + as qjj,{l - 75)?. u^l - 75)c] , 

(15) 

accompanied by the emission of the virtual photon. Here qij denote the d or s quark fields. 
The coefficients ai = 1.2 and 02 = —0.5 have been determined from the experimental data 
on nonleptonic charm meson decays in the extensive analysis based on the factorization 
approximation of p^. We also systematically undertake the factorization approximation to 



evaluate the matrix element for the product of the currents (15). 

In order to treat the transition among physical particles, we shall use the effective La- 
grangian approach with heavy pseudoscalar D, heavy vector D*, light pseudoscalar P and 
including also light vector V degrees of freedom. The later are necessary, since they play a 
dynamical role in the photon emission from a meson via vector meson dominance (VMD) 
and lead to the resonant spectrum in terms of invariant di-lepton mass mu. We organize 
various effective interactions among the mesonic degrees of freedom following the Heavy Me- 



son Chiral Lagrangian approach p8[, which is reviewed in and is most likely the best 
suited framework for treating the problem under investigation. It embodies two important 
global symmetries of QCD: the heavy quark spin and flavour symmetry SU{2Nf) in the 
limit mc 00 and chiral symmetry SU{3)l x SU{3)r, spontaneously broken to SU{3)v, in 
the limit mu^d,s 0. The light vector mesons are introduced by promoting the symmetry 

G = [SUi3)L X SUi3)R],lobal/[SUi3)v]global to G' = [SU{3)l X SUi3)R]global X [SUi3)v]local, 

where the light vector resonances are identified with the gauge bosons of [SU{3)v]iocai [pOj - 
One is free to fix the gauge of [SU{3)v]iocai and the two theories, based on the groups G and 



G', are equivalent up to terms with derivatives on the light vector fields |^ 



Keeping only the kinetic and interaction terms of the lowest non-trivial order, the La- 
grangian has the form pU|, 



C = -4|tr[^,^1 + a tr[iV, - p,)^] }+ J^tr[F,.(p)F^'^(p)] 

+ tTr[H,v^{Si,ad^ - ileSbaA^' + V^, - «:(V^ - p'')6a}i?a] + igTr[H,^^-f,AlH,] , (16) 



with 



A^ = lieid, + teQA^)^ - + teQA^)^] 
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Q = diag(2/3, —1/3, —1/3) and photon field A^. The hght fields are incorporated in 



/4 



^ = exp 



f 
I 



V2 



\ 



K- 



71 



jlL _i_ \2R. _i_ ^0 1 



K+ 

2->?8 



±28 _|_ i2o.l , 



where 773 and r^o contribute to r] — rj' mixing as in |]T3|] with = —20 ± 5°. 
pseudoscalar and vector fields of fiavour cqa are incorporated in 



Ha = I'H^l' . 



(17) 



The heavy 



Above, / = 132 MeV is the pseudoscalar decay constant and gy = 5.8 is the VPP 
coupling p9| , pO| . We fix a = 2 assuming the exact vector meson dominance, when the light 



pseudoscalars interact with the photon only through the vector mesons p9| , pO| , p2| . We shall 
use g = 0.59 ± 0.06, obtained by CLEO from the measurement of the widths D*^ D^tt'^ 
and D*^ D^n^ [^. The parameter n will eventually turn out to be multiplied by a small 
factor mp in the D —>■ Pl^l^ amplitudes and its contribution is negligible. 

The bosonized weak current coming from the light quarks is obtained by gauging (|T^) 



ba 



(19) 



The weak current ga7'^(l — 75)c transforms under chiral SU{3)l x SU{3)ji transformation 
as (3^,, In) and it is linear in the heavy meson fields and Z^*" [^, |3^ 

qa^a - 7')c ^ I^/dv^ Trira - 75)^^6^11 (20) 
+ a^Tr[r,H,{p^ - V^bcC] + a^Tri^^.H.v^ip^ - V"),,eL] + - 

This current is the most general one at the leading order in the heavy quark and next- 
to-leading order in the chiral expansion. The parameters ai and ^2 are determined from 
experimental data on Br, T i/Tt and r+/r_ of the decay K*^e^V(. ||13[|. Among the 

eight sets of solutions for three parameters [^], we use the set ai = 0.14 ± 0.01 GeV^/^ and 
a2 = 0.10 ± 0.03 GeV^/^ which agrees with the measured form factors. 

We shall calculate a larger group of -D — Pl^l~ decays, rather than only those related 
to c — i> ul^l~ transition. The list of decays considered is given in Table |^. The Feynman 
diagrams for the long distance contributions to D — > Pl^l~ within our framework are given 
in Fig. H The Lagrangian (|15D contains a product of two left handed quark currents, each 
denoted by a dot in a box. We organize different diagrams according to the factorization of 
the non-leptonic effective Lagrangian ([T5|): 



The long distance penguin contribution in Fig. Qa is induced by [57^5 — 
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/+ /+ /+ /+ 




p p p p 

(a) Long distance penguin contribution. 

/+ /+ /+ /+ /+ Z+ 1+ 




p 

(b) Long distance weak annihilation contribution 



Figure 2: Long distance contributions to D ^ Pl^l~ decays. The vector meson V denotes 
p°, uj or 0. The box denotes the action of the nonleptonic effective Lagrangian ([T5|). The 
box contains two dots each denoting a weak current in the Lagrangian (^). 

• The long distance weak annihilation in Fig. is induced by a product of the 
weak currents, where one current has the flavour of the initial D meson, while the other 
has the flavour of the final P meson. Vector resonances do not enter as intermediate 
states R in the weak transition D ^ R followed hj R —>■ P7* or D — > Rj* followed by 
weak transition R ^ P, since parity is conserved in D P'-f* process. 

The Lagrangian (|TBp and the weak currents (JT^, ( PD| ) are invariant under the electro- 
magnetic gauge transformation and automatically lead to the gauge invariant amplitude 
of the form (|13D . This is due to the fact that the vector field and the vector current 

= ieQA^ + I (^^9^^ + ^<9^^^) always appear in the gauge invariant combination V^^ — Pfj, 
and the resonant and nonresonant diagrams in Fig. ^ come in pairs. 

We incorporate SU{3) symmetry breaking by using the physical masses, widths and decay 
constants, given in Tables ^ and § of Appendix with the following definition 

(0|j^|P) = tfpp^, {0\f\D) = -if op'. (0|/|V) = gve', {^\f\D*) = ifo^mo^e' (21) 

and properly normalized = gi7'^(l — 75)5'2- The assumptions for extrapolating the am- 
plitudes away from where the chiral and heavy quark symmetries are good, are discussed in 
Appendix. The amplitudes for the diagrams in Fig. ^ are given by Eq. (|26|). 

4 The results 

The allowed kinematical region for the di-lepton mass mu in the D Pl^l^ decay is 
mil = ['2''mi,m£,—mp]. The long distance contribution has resonant shape with poles at mu = 
mpO,m^,m^. There is no pole at mu = since the decay D P'-f is forbidden. The short 
distance contribution is rather flat. The spectra of D Pe^e~ and D — > Pp^ p~ decays 
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D - 
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I = H, e 


BrsM ^ Br'^'' 
I = fi, e 


I = e 




I = fi 




D'- 


^ K^l+l~ 





4.3 X 10-^ 


< 1.1 X 10" 


-4 


< 2.6 X 10- 


-4 


Dt- 


n+l+l- 





6.1 X 10-6 


< 2.7 X 10- 


-4 


< 1.4 X 10- 


-4 


- 


TiH+l- 


1.9 X 10-9 


2.1 X 10-^ 


< 4.5 X 10- 


-5 


< 1.8 X 10- 


-4 


- 


r]l+l- 


2.5 X 10-1° 


4.9 X 10-s 


< 1.1 X 10- 


-4 


< 5.3 X 10- 


-4 


- 


r]'l+l- 


9.7 X 10-12 


2.4 X 10-1° 


< 1.1 X 10- 


-4 


< 5.3 X 10- 


-4 


D+ - 


Tc+l+l- 


9.4 X 10-9 


1.0 X 10-*^ 


< 5.2 X 10- 


-5 


< 7.8 X 10- 


-6 


Dt- 


^ K+l+l- 


9.0 X 10-1° 


4.3 X 10-* 


< 1.6 X 10- 


-3 


< 1.4 X 10- 


-4 


D+ - 


^ K+l+l- 





7.1 X 10-9 


< 2.0 X 10- 


-4 


< 8.1 X 10- 




D'>- 


^ KH+l- 





1.1 X 10-9 











Table 2: The branching ratios for nine D Pl^l^ decays in the standard model. The short 
distance contributions, induced by the c —>■ ul^l~ transition, are given in column 2 and are 
small. The total branching ratio is therefore dominated by the long distance contribution 
and is given in column 3. The experimental upper bounds are given in the last two columns 
p!3| , |Tl], |I^: the E791 analysis |]1T| considers and decays, while the new analysis of 
FOCUS ||12| considers only decays. The MSSM has insignificant effect on the total rates 
of D Pl+l- decays. 



in terms of mu are practically identical. The difference in their rates due to the kinematical 
region mu = [2me, 2m^] is small and we do not consider them separately. The predicted 
branching ratios for nine decays in the standard model are given in Table ^ together with 
the available experimental data ||TT| , |l^, |l3l. The short distance contribution, as predicted 



by the standard model, is given in the second column and is small. The total branching ratio 
is therefore dominated by the long distance contribution and is given in column 3. 

The differential branching ratio dBr /dm\ for the Cabibbo allowed decay n^l^l^ , 

which arises only via the weak annihilation, is presented in Fig. |^a. In Fig. ^3, we present the 
Cabibbo suppressed decay D 7il~^l~ , in which the kinematical upper bound on di-lepton 
mass mj^"^ = mD — mp is the highest. The dashed and dot-dashed lines denote the long and 
short distance parts of the rate in SM, respectively, while the solid lines denote the total rate. 
The long distance contribution decreases in the kinematical region above the resonance and 
the short distance contribution becomes dominant. Thus, the decays D~^'^ 7?+'°/"'"/- at 
high mil might present a unique opportunity to probe the flavour changing neutral transition 
c ul^l~ in the future. As the pion is the lightest hadron state, this interesting kinematical 
region is not present in other D — * Xl^l~ decays. 

The differential distribution for ti^I^I~ , given in Fig. |^, indicates that the high 

di-lepton mass region might give an opportunity for detecting c — > ul^l~ . Before making 
a definite statement on such possibility, we should examine this kinematical region of high 
di-lepton mass m D nl^l^ decays more closely. For instance, in this region the excited 
states of the vector mesons p, u and may become important. We attempt a rough estimate 
of the additional long distance contribution arising from the first radial excited states pi, 
Ui and 01 (^5"!) and first orbital excited states p2, oji and 02 (^-Di). The knowledge of 
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0.0 0.5 1.0 1.5 2.0 2.6 3.0 3.5 0.0 1.0 2.0 3.0 

m^JGeV"] m'^JGeV'l 

(a) (b) 

Figure 3: The differential branching ratios dBr /dm\ as a function of the invariant di-lepton 
mass mfi for the Cabibbo allowed decay — >■ vr+Z+Z" (a) and Cabibbo suppressed decay 
—¥ 7r+/+/^ (b). The dashed line denotes the long distance contribution, the dot-dashed 
line denotes the c ul^l~ induced short distance contribution, while the solid line denote 
the total standard model prediction. The 7i^l^l~ arises only via the long distance 

contribution. 



their masses, decay widths and couplings to other particles is poor at present. We use the 
measured masses and widths, taken from [|1^, ^] and compiled in Table ^. Due to the lack 
of the experimental data on the leptonic decay widths we use the magnitudes of the 



decay constants gy as predicted by the quark model in ||36|lP| and compiled in Table |^. At the 
same time, we assume that the excited vector mesons couple to the charmed mesons with 
the same couplings as the corresponding ground state vector mesons p, u and (p. In this case, 
the corresponding amplitudes (pB]) are obtained by replacing the coefficients A^^i and Mi by 
the expressions given in (^81) . The differential branching ratios for D —>■ 7r/i+/x~ decays are 
given in Fig. ^. The thick and thin dashed lines denote the long distance contributions with 
and without excited vector mesons, respectively. The short distance contribution, denoted 
by dot-dashed line, is still dominant in the kinematical region of high mu in spite of the 
excited vector resonances. 

The possible enhancement within the general MSSM, discussed in section 2.1, is presented 
in Fig. ^ and is probably too small to be observed in any D Pl^l~ decay. The solid lines 
represent the standard model prediction for the D — > 7rl~^l~ branching ratios. The dot- 
dashed lines represent the best enhancement in the general MSSM and indicate that the 
D Pl^l^ rates are rather insensitive to the large supersymmetric enhancement of cj. The 
value of cj is manifested in c ^ ul^l~ at small mu (see Eq. and Fig. |1]), while it's effect 
is suppressed in D — >■ Pl^l~ decays do to the factor in the general expression for the 
D FY amplitude (|T3p. 



"^The decay constant /y, defined in is related to gv, defined in (|2]]), by: /p V^iripfp, f^; 
SV^nii^fu, and -3m^f^. 
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1.5 



Figure 4: The differential branching ratio for D — > Tr/i"*"//" decays. The thick dashed hnes 
present the long distance contribution incorporating the ground state and the excited vector 
mesons. The thin dashed lines present the long distance contributions due only to the ground 
vector mesons. The short distance contribution, denoted by dot-dashed line, is dominant in 
the kinematical region of high ma, in spite of the excited vector resonances. 




Figure 5: The biggest possible enhancement of D — *• Tr/i+zi" rates within the general MSSM, 
discussed in section 2.1, is denoted by the dot-dashed lines. The solid lines represent the 
standard model predictions. The effect of supersymmetry is screened by the uncertainties 
present in the determination of the long distance contributions and is probably too small to 
be observed. 



5 Conclusions 

We have presented the first predictions for rare charm meson decays D — >■ Pl^l~ with 
P = 7r,K, 7], 7]' in all nine possible channels; the previous analysis has considered only 
the D — >■ ttI^I' channel. The long distance contributions are found to dominate over the 
short distance contributions, which are induced by c — > ul^l~ in the Cabibbo-suppressed 
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decays. We have used the theoretical framework of Heavy Meson Chiral Lagrangian with 
the recently determined value of the strong coupling g from the measurement of D* Dtt 
width. Our predictions are compiled in Table 0. The decay tt^I^I~ is predicted at 

the highest branching ratio of 6 x 10~^. The best chances of the experimental discovery are 
expected for —>■ n^l^l^ , which is predicted at 1 x 10~^ and has the upper bound 8 x 10~^ 
at present. The limits on and modes at the level 10~^ are expected from CLEO-c 



and B-factories, while the limits on Dj^ modes are expected to be an order of magnitude 
milder 

The only possibility to look for c — ul^l^ transition is represented hj D -kI^I^ decays 
in the kinematical region of mu above the resonance 4>, where the long distance contribution 
is reduced (see Fig. ^). 

We have explored the sensitivity of the c ul^l~ within two scenarios of physics beyond 
SM. The effect due the exchange of the flavour changing Higgs in Two Higgs Doublet model is 
found to be negligible. The general Minimal Supersymmetric Standard model can enhance 
the c — > rate by up to a factor of three (see Table 0). This effect is due to the 

large supersymmetric enhancement of cy and is seizable at small mu in c ul^l~, but it 
is unfortunately very small in the hadronic process D Pl^l~ as the decay D P7 is 
forbidden (see Fig. 

The kinematics of the processes D — > Vl^l^ would be more favorable to probe the possible 
supersymmetric enhancement at small m//, but the long distance contributions in these 
channels are even more disturbing |0. The large supersymmetric enhancement of the Willson 
coefficient Cy is manifested in c — * ^7 decay and can enhance the standard model rate ~ 10~^ 
by up to two orders of magnitudes 0]. Such enhancement could be probed by observation 
of B(. — s> 5*7 or by measuring the relative difference Br{D^ p^'-f) — Br{D^ owf) [Q. 

6 Appendix 

The short distance part of the D Pl^l^ amplitude, induced by the transition c — > ul^l~ , 
contains the following form factors 

(P(p')|g7^(l-75)c|/^(p)) = ip + p')Mq') + ip-p')J-iq') 
{P{p')\qa^u{l ± l5)c\Dip)) = is{q^) [{p + p')^q, - q,{p + p% ± te^^x^P + p')V] (22) 
defined using operators in (|I]). The short distance amplitude is then given by 



^SD[B(p) ^ p(p_ ,)(+(-] = , S^eX',Vu. 



u{p^)]iv{p+) , 
(23) 

where we neglected the nearly vanishing cio, Cg and coefficients in SM (^) and MSSM (|^). 
In the heavy quark limit, the form factor s can be expressed in terms of the form factors f± 
at zero recoil ||39(] p| and we assume the relation to be valid for all q^ 

^(?') = ^[/+(?') ■ (24) 



^This relation is not written correctly in l39| and is corrected in [E9l 
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The semileptonic form factors f± in the Heavy Meson Chiral Lagrangian approach, extended 
by assuming the polar shape, are given by |]28| , |29| p| 



-K 



D 



DP 



9 



rriD — Trip 



m 



D' 



^max 



m 



(25) 



D'* 



with Kdp given in Table 

The long distance amplitude is given by the diagrams in Fig. ^ The long distance 
penguin diagrams in Fig. |]a are expressed in terms of the form factor (^Sf). The weak 
annihilation contribution in Fig. ^ is determined by assuming that the vertices do not 
change significantly away from the kinematical region, where the heavy quark and chiral 
symmetries are good. We expect this to be a reasonable approximation in D meson decays. 



At the same time we use the full heavy meson propagators 1 / (pj^ ■ 



■m 



instead of the HQET 



propagators l/{2mvk) In the limit nip <^ m^i, the bremsstrahlung-like diagrams in Fig. 

cancel exactly, as explained in detail in the Sections 3.3.3 and 5.5.1 of Only the non- 
bremsstrahlung weak annihilation diagrams in Fig. ^ja render the non- vanishing contribution. 
The long distance amplitude is given by ||32[ 



\LD 



(26) 



annih. \\l J \^ ' ^ annih. \\l J 

hremsstrahlung non — brem. 



annih. \\l 



peng 

annih. 
hremsstrahlung 

\LD ! „2 

annih. 
non — brem. 



{q') = f!^Li,Mr>{qV. 



m. 



mj^ — nil 



- vVolai ^ a2) 



9v_ 
V2 



with Cabibbo factors /(^^^f,^ and the coefficients Mi(g^) and K^j^p as given in Table ^ The 
coefficient Ni equals 



Ni{q' 



9l 



-ml + iVpnip ?,{q'^ - + iV ^m^) Slq"^ - ml + iT^m^) 



ml 3m,^, 3m 



2 ' 



while the coefficients m[^^ are given in terms of M f , Mf and M^^ in Table | 



DO 



9p 



9p 



— rn? + iV 3(g^ — + zF^m^) ^^rn'^ 



D+ 



9p 



+ 



9p ^ 9u. 



g2 _ ^2 _|_ iYpirip 3(g^ — + iT^m^) m^ 3m| 



Ml 



D^, 



3(g2 -ml + iV^m^) 3m 



2 ■ 



(27) 



^ Different form factors f± were used together with g ~ 0.27 in |3|]. These form factors would overproduce 
the semileptonic decay rates for the value g ~ 0.59 recently measured by CLEO [p3[. 
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Note that A^i(O) = Mi(0) = for r(0) = and there is no pole arising from to the photon 
propagator at = 0. The relative sign of the short and long distance penguin amplitudes 
agrees with the Ref. which is based on assumption of quark-hadron duality. 



i 


D Pl+l- 


JCabb 




^DP 


1 











2 


7~\ J- 1 T 1 T 









3 




-a2VudV:a 




1 


4 


7]l+l- 




-^fMf cos^p 


COS Op sin Bp 
V6f V3f 


5 


7]'l+l- 




-y/fMf° sin^p 


sin 9 p 1 cos 9 p 

Vef ' V3f 


6 


D+ ^ n+l+l- 


—aiVudV*^ 




1 


7 


D+ K+l+l- 


aiVudV:^ 




1 

/a- 


8 


D+ K+l+l- 




Mr 





9 


^ K^l+l- 


-a2Vusv:, 








Table 3: The Cabibbo factors fQ^M, the coefficients Kj^'p and the functions for nine 
D Pl+l- amphtudes in (|2|). 



H 


niH 


fH 


P 


nip 


fp 




[GeV] 


[GeV] 




[GeV] 


[GeV] 


D 


1.87 


0.21 


TT 


0.14 


0.135 


Ds 


1.97 


0.24 


K 


0.50 


0.16 


D* 


2.01 


0.21 


V 


0.55 


0.13 








v' 


0.96 


0.11 



Table 4: The values of the meson masses, decay constants and decay widths [0. The 
measured decay constants fo and fw have sizable uncertainties and the values are taken 
from lattice QCD results [|3^ . 



In order to account for the contributions of the excited vector mesons pi,2, ^1,2 and 01,2, 
as described in the main text, the coefficients A^i and Mi are replaced in the formulas above 
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P UJ (f) 


Pi t^i 01 


P2 ^2 02 


m[GeV] 
r[GeV] 
9v[GeV^] 


0.77 0.78 1.0 
0.15 0.0084 0.0044 
0.17 0.15 0.24 


1.45 1.46 1.69 
0.31 0.24 0.3 
0.11 0.11 0.23 


1.66 1.66 1.88 
0.4 0.1 0.3 
0.07 0.07 0.12 



Table 5: The masses, widths and decay constants of ground |]T3| and excited |^ ^ vector 
mesons. 



^ by 



A^i ^ A^i + J] 



9lu 



+ 



k 

3pk 9uj^ 



Hk 



m 



3m2,. 3m^^ ' 



Pk 

2 



Mi 



DO 



9pk 



+ 



9^k 



9pk _|_ 9ujk 



k=l 
2 



D+ 



9pk 



+ 



9u!k 



9pk _|_ 9ujk 



^^q^ -ml^+iTp^mp^ 3{q^ - ml^ + iV^^m^J m\ 3m; 



2 ' 

'^k 



'^94>k 



1 3(g2-m2^ +zr<^,m^J 3m; 



2 

4>k 



(28) 
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